Charge 2e skyrmions in bilayer graphene 
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Quantum Hall states that result from interaction induced lifting of the eight-fold degeneracy of 
the zeroth Landau level in bilayer graphene are considered. We show that at even filling factors 
electric charge is injected into the system in the form of charge 2e skyrmions. This is a rare example 
of binding of charges in a system with purely repulsive interactions. We calculate the skyrmion 
energy and size as a function of the Zeeman interaction, and discuss signatures of the charge 2e 
skyrmions in the scanning tunneling microscopy experiments. 



Introduction. The four-fold valley and spin degeneracy 
of Landau levels (LL) in monolayer and bilayer graphene, 
the recently discovered two-dimensional semimetals [l], 
HI, 13 1 gives rise to interesting phenomena at high mag- 
netic fields, where the Coulomb interactions between 
the electrons become important. In the monolayer, the 
Coulomb interactions lift the LL degeneracy, giving rise 
to new spin and/or valley polarized incompressible quan- 
tum Hall (QH) states [!, 0, [|| . The Hamiltonian of the 
interaction induced quantum Hall states is approximately 
SU(4) symmetric UQ with respect to the rotations in 
the combined spin/valley space. The splitting of the LLs 
thus corresponds to the spontaneous symmetry break- 
ing of the 5'J7(4)-symmetric quantum Hall ferromagnet 
(QHFM). The precise order in which spin and valley de- 
generacy get lifted is determined by the interplay between 
the Zeeman interaction and valley anisotropy [^, Hcj |. 
both of which are much smaller than the Coulomb in- 
teraction. The spin- and valley-polarized QH states were 
predicted 11| to feature spin and valley skyrmions, which 
are smooth topologically nontrivial textures of the ferro- 
magnetic order parameter that carry the electron charge 



12j. The QHFM states in the monolayer and bilayer 



graphene are also expected to have interesting edge states 
properties 13, 14, 15 1, as well as unusual spectrum of the 
low- lying collective excitations 11, [l(|. Also, textured 
states in the vicinity of integer filling factors have been 
predicted [llj ]. 

Bilayer graphene features a LL at zero energy, which 
has a twofold orbital degeneracy: in each valley there are 
two zero-energy states (a = 0,1), with wave functions 
corresponding to the ground state and the first excited 
state of the magnetic oscillator [l7| ■ Taking into account 
valley and spin degeneracies, the zeroth LL in the bi- 
layer is eight-fold degenerate. Coulomb interactions are 
expected to lift the eight-fold degeneracy [3]. In this 
paper, we consider the interaction induced QH states at 
even filling factors, and analyze their new properties aris- 
ing due to the orbital isospin. We shall sec that these QH 
states exhibit interesting collective and topological exci- 
tations. We predict that pairs of the excitations of charge 
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FIG. 1: (a) Bilayer graphene lattice. Perpendicular elec- 
tric field E generates effective valley Zeeman interaction, 
A„ = eEd, where d — 0.34 nm is the separation between the 
layers, (b) The order of the zeroth LL splitting, assuming that 
effective valley Zeeman interaction A„ favors K valley states, 
and A v < E z . (c) Texture corresponding to the charge 2e 
skyrmion at v = ±2. Vectors illustrate the rotation of the 
order parameter in the valley space. 



e bind into skyrmions that carry charge 2e. Such bind- 
ing of charges is surprising, because the Coulomb inter- 
actions between electrons are purely repulsive. Another 
example of such binding was predicted to occur in the 



spin QHFM with small Zeeman interaction [18|, 1191. The 
weak pairing of skyrmions considered in Refs. [La, [l9| . 
however, can occur only when the Zeeman interaction is 
extremely small; in contrast, charge 2e skyrmions in bi- 
layer graphene can be thought of as robust tightly bound 
pairs, which exist in a wide range of the effective Zeeman 
interaction. 

Below we analyze the dependence of skyrmion energy 
and size on the potential difference between the two lay- 
ers of bilayer graphene, A„ = eEd (see Fig. la), which 
favors one of the valleys and therefore acts as a valley 
Zeeman interaction. Owing to the fact that the surface 
of bilayer graphene is exposed, the skyrmion size can be 
measured using scanning tunneling microscopy (STM). 
Furthermore, we find that slightly away from even filling 
factors, |Ai/| = \v — 2M\ <C 1, there is a finite density of 
charge-2e skyrmions in the system. At small density, the 
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skyrmions form a triangular lattice, while above a critical 
density, , they form a bipartite square lattice 23], 24 1 . 
The phase transition between the two lattice structures 
can be observed by STM. 

The effective Coulomb interaction Hamiltonian for the 
zcroth LL in the bilaycr is approximately SU (4) sym- 
metric in the valley-spin space, however, the symmetry 
in the orbital isospin space is broken due to the differ- 
ent orbital wave functions of the two states [lij ]. This 
results in the following picture of the zeroth LL splitting: 
at even filling factors (v = 2M filled sub-levels) M pairs 
of orbital states with the same valley and spin are filled, 
while the states at odd filling factors, v = 2M + 1 are 
obtained from the v = 2M QH state by filling one of the 
remaining states with orbital isospin a — 0. This order 
of the zeroth LL splitting is due to two facts: (i) ex- 
change energy within the LL with isospin a = is higher 
than that for the LL with isospin a = 1; (ii) there is ex- 
change energy between filled a = and a = 1 LLs with 
the same spin and valley, which makes the energy of the 
state where a = 0, 1 LLs with the same spin and valley 
are filled (e.g. OK | and IK f ) lower than the energy of a 
state polarized in the orbital space along a = direction 
(e.g., OK |, OK' "f LLs are filled). The order in which 
valley and spin degeneracies get lifted is determined by 
the competition between the symmetry-breaking terms: 
the Zeeman interaction, E z — g/isB, and the effective 
valley Zeeman interaction A,.. In the experiment A,, is 
typically small, and it can be tuned by gates [22j. We 
assume that A„ is tuned to be smaller than E z . Fur- 
thermore, for simplicity, we assume that A„ is small but 
non-zero and favors the K valley 25]. This leads to the 
splitting picture illustrated in Fig. lb. In the following, 
we shall be especially interested in the states at filling 
factors v = —2, +2, marked by arrows in Fig. lb. Since 
these two states are related by the particle-hole symme- 
try, we shall focus on the state at v = — 2. 

Landau levels in bilayer graphene. We start with 
recalling the Landau level spectrum in the bilayer 
graphene [rH |. The low-energy excitations near the 
K, K' point are described by the Schroedinger equation 
etpK,K' = HK,K'tpK,K' , with the Hamiltonian given by 
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where the effective mass m can be expressed in terms 
of the interlayer hopping amplitude 71 0.39 eV and 
the Fermi velocity in the monolayer vp 1=3 10 6 m/s, 
m = 7i/2-y|,. For the K valley the upper (lower) compo- 
nent of the wave function corresponds to the amplitudes 
on the sublattices A(B) (see Fig. la), which belong to 
different layers. For the K' valley the order of compo- 
nents is reversed, such that the upper (lower) component 
corresponds to the B{A) amplitude. 

To analyze the LL spectrum, we choose the Landau 



be classified according to the value of the wave vector k y , 
*l J K,K'(x,y) = e lkyV ipK.K'(x). The wave vector k y trans- 
lates into the guiding center position, X — k y £ 2 B , where 
£b = \JficjeB is the magnetic length. Below for simplic- 
ity we shall choose units where £b = 1- The effective Id 
Hamiltonian for ipK,K' (%) takes the following form, 



H k ,k> = -tu^c 
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ax = i{d x + (x- X)), 



m 

where lo c = eB/mc is the cyclotron energy. The Hamilto- 
nian ^ has two zero modes with the following wave func- 
tions, ^ K . K >(x,y) = e tXy (0,Lp a: x(x)) 7 a = 0, 1, where 
(f a x{x) denotes the a-th excited level of the magnetic os- 
cillator. Below we shall denote the annihilation operators 
of the zero modes by c a K jf , k — K, K' . 

Coulomb interaction. Now we proceed to the analysis 
of the zeroth LL splitting. We neglect the LL mixin g (th e 
effects of LL mixing will be considered elsewhere [23|L 
which allows us to project the Coulomb Hamiltonian onto 
the zeroth LL. The interaction Hamiltonian can be writ- 
ten in the following form, 



(3) 



L x L y is the 



where p(q) are the density operators, S 
sample volume, the valley indices, and the ma- 

trix element is given by V(q) = ^p- (26| . 

In order to project the Coulomb interaction onto the 
zeroth LL, we introduce projected density components as 
follows, 



PM) = Y] F ab (<lK\ pf (q) = V exp (zq x X 
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were X± = X ± f , and F 00 (q) = e~« '\ ^(q) = 
(1 — q 2 /2)e~ q / 4 are the usual form-factors for the lowest 
LL and the first excited LL, and F i(q) = -V«±^ e -<i 2 : 

Fio(q) = 9v ft* e~ q / 4 , are the form-factors correspond- 



ing to the density components which mix the two orbital 
states. The effective Coulomb interaction within the ze- 
roth LL is obtained by plugging Eq.(j4]) into Eq.([3|). 

Nature of the states with even filling factors. We now 
analyze the order in which the eight-fold degeneracy of 
the LL gets lifted. The split QH states with filling factor 
\v\ < 3 (^+4 filled sub-levels) correspond to the following 
wave functions, 
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gauge, 



Bx, A x = 0, for which the eigenstates can 



where S are linear combinations of the operators: 

Ax = E uUA, K , s ,x, (6) 
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where U is a unitary matrix, and s is the electron spin. 

To find the ground state for i/ = -2, we compare the 
energies of two states: (i) two a = LLs with different 
valley and/or spin indices are filled, for example, d\ = 
c o k t ' ^2 = c o K' t ' an< ^ ( u ) a — an d a = 1 LLs with 
the same valley and spin indices are filled, d\ — K 

c?2 = cj ^ | . The energy of the first state is twice the 
exchange energy of a non-degenerate lowest LL, 

(^ nt)l = -2iVA 0) *o = \yJ\^ (7) 

where TV is the total number of states in one non- 
degenerate LL. Averaging the effective Coulomb inter- 
action over the second state, we obtain 

(H int ) 2 = ~NA . (8) 

Thus the energy of the second state is lower than the 
energy ([7]) of the first state, and the spin- and valley- 
polarized state |-0o) is the ground state at v = —2. The 
state at v = +2 can be obtained from \if; ) by charge 
conjugation. 

Charge 2e skyrmions. Now we proceed to discussing 
excitations of the v — —2 QH state. The lowest-energy 
electron- hole pair at v = — 2 is obtained by removing an 
electron with orbital isospin a = 1 from the filled LL, 
and putting it into one of the empty LLs. The energy 
of such excitation, E e h = |Ao, is lower than the energy 
E' eh ~ 4Ao of a particle-hole excitation that is obtained 
by removing an electron with isospin a = 0. 

In some QHFMs, the lowest energy charge excitations 
are skyrmions, which are topologically nontrivial smooth 
textures of the order parameter [l2|. On the qualitative 
level, the textures carry charge because the charge and 
spin (and /or valley) dynamics in the QHFM are entan- 
gled EiO. 

Can skyrmions exist in bilayer graphene? Skyrmions 
of charge e are energetically unfavorable because they 
involve flipping valley isospin (or spin) for either a = or 
a = 1 states in some region, and in that region the filled 
a = and a = 1 states would have different valley isospin 
(spin) , which leads to a loss of the exchange energy Ao 
per flipped valley isospin. 

Another possibility is skyrmions of charge 2e, which 
can be created by making two identical valley textures 
for a — and a — 1 orbital states. Such textures are de- 
scribed by a unit vector n, with n z = — 1(+1) correspond- 
ing to filling K(K') states. On the intuitive level, we 
expect such textures to be energetically favorable: since 
a = and a = 1 orbital states rotate simultaneously, no 
exchange between and 1 states is lost. Below we find 
the energy of the 2e skyrmion, and, by comparing it to 



the energies of the single-particle excitations, establish 
that such skyrmions are indeed energetically favorable. 

Before we proceed to the quantitative analysis of 
charge 2e skyrmions, we would like to compare excita- 
tions at the even and odd filling factors. For simplic- 
ity, let us consider the excitations of the state v = —3, 
which corresponds to filled OA | LL. The lowest energy 
electron-hole pair is obtained by removing an electron 
from the OA | LL and putting it in the IK f LL; ow- 
ing to the exchange between OA" f and 1A" | states, the 
energy of such excitation, E°£ d = Ao, is lower than the 
energy E°£ d = 2Aq of an excitation where the excited 
electron resides in a LL with a different valley/and or 
spin index. The existence of orbital skyrmions at v = — 3 
is unlikely, because such skyrmions correspond to filling 
a = 1 states in some region, which leads to a loss of 
the exchange energy equal to Ao/4 per flipped orbital 
isospin. 

Skyrmion energy. In order to compute the energy of 
charge 2e skyrmion, we derive an effective Hamiltonian 
describing textures of the order parameter, 

\^)=e- id \^ Q ). (9) 

In our analysis, we follow the microscopic approach devel- 
oped in Ref. [2fJ; as we shall see below, the dynamics of 
the order parameter in the bilayer graphene is richer than 
that in the case of SU(2) and Sf7(4)-symmetric QHFM, 
owing to the presence of the orbital degree of freedom. 
We parametrize the rotation operator O as follows, 

0= £ ^(q)^(-q), (10) 

^(-q) =T, elQ ^ T -f Lc L,x^,x + , (11) 
x 

where r are the Pauli matrices. The rotation (fTU)) is 
described by four complex parameters, 

Ua = Ka+i^la, O = 0, 1, « = ^ +iQ 10 , W = ^l+^Ol' 

(12) 

The parameters u (ui) correspond to rotations that in- 
volve OA and OA' (1A and 1A') states, while v, W 
parametrize rotations which transform OA into 1A', and 
1A into OA', and vice versa. To simplify calculations, 
we assume that the rotations are small (\u a \ <§C 1, |u| <§; 
1, |w| <^ 1). Then we can expand the texture energy 

E= {4>o\e l ° He~ l °\ipo} — (tpo\H\ip ) in series in the pow- 
ers of O. Doing that and evaluating the mass and stiffness 
terms, we obtain (23l |. 
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£ , = / _L£[40 u *0 Uo + 7fl«;Sui + 6{5uS9wi + c.c.} + 2\/2{uS(5w + dw) + c.c.} + \/2{^(<9u + dw) + c.c.} 



Po 
2 



2/ 



-2{(%*<9i; + <9w*<9u>} — 3{(%*<9w + c.c.} + 2(itg — u*)( u o — «i) + 4v*v + 3w*w], p 



(13) 



r 



We are interested in the low-energy excitations, where 
a = and a — 1 states are rotated simultaneously 
in the valley space (see Fig.lc). This corresponds to 
setting Mo = u\ = u, which ensures that the mass 
term that contains uq,ui components in Eq. (|13p van- 
ishes. Minimizing the energy with respect to v,w yields 



3\/2 



du 



= V2di 



where ellipsis denotes 



the higher order gradient terms. Notice that v, w are 
proportional to the gradients of it, which implies that 
v, w are much smaller than u for the case of slowly vary- 
ing textures. Substituting the expression for v,w into 
Eq. (|13[) . we obtain the gradient term for it, 



E s t — —2p s 



dzdz- 25 

du du, p s = — po- (14) 



For what follows, it is convenient to rewrite the stiffness 
energy (fl"4|) in terms of the 0(3) order parameter n = 
(-u y ,u x ,0), 



E st = ^f d 2 r{d,nf 



(15) 



Although we have derived the above equation assuming 
that n deviates slightly from n = (0,0,-1), due to the 
rotational invariance in the valley space Eq. p5p is valid 
for an arbitrary slowly varying configuration of the order 
parameter. 

As our next step, we evaluate the charge density of the 

texture, Sp = (V>o|e l0t /3e~ l0 |?/;o) — (^oIpI^o), where p is 
the density operator. We find that the charge density is 
twice the Pontryagin index density, 

5p(r) = 2ep(r), p(r) = --i-e^n^n x cVn]). (16) 

Notice that this relation differs from the usual SU(2) 
QHFM case 12, 2^] by a factor of 2, which corresponds 
to the fact that the texture rotates states in both a = 0, 1 
LLs. 

Apart from the stiffness term (|15[) , there are two other 
contributions to the texture energy: the valley Zeeman 
term and the long-range Coulomb interaction, 



H z = A v n J d 2 rn z , (17) 

H coul = \ [ d 2 vd 2 v'8p{r)—^—5p{v'), (18) 
2 J |r-r'| 

where uq — l/2ir£ 2 3 is the LL density of states. 



The simplest topologically nontrivial texture of the or- 
der parameter n has topological charge 1 and an elec- 
tric charge ±2e. This is to be contrasted with the 
usual skyrmions [l2j], which carry charge ±e. In the 
limit of vanishing A„ , the Coulomb repulsion (fT5)) forces 
skyrmions to be infinitely large, l s — > oo, where l s is the 
skyrmion size. Then the skyrmion energy is determined 
solely by the stiffness term, 



25 

E sk = 4np s = — A . 
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(19) 



The energy of the skyrmion-antiskyrmion pair, 2E s k = 
25Ao/8, is lower than the energy of two electron-hole 
pairs, which equals 7Ao. Therefore, in the limit A„ — > 0, 
pairs of electron (hole) excitations bind into charge 2e 
skyrmions (antiskyrmions) . 

At finite A^ the skyrmion size is determined by the 
competition between the effective valley Zeeman and 
Coulomb energies Optimizing the skyrmion energy 
with respect to its size, we find with logarithmic preci- 
sion, 



(20) 



where A^ = A w /(e 2 /£s). The skyrmion energy is in- 
creased compared to the case A„ = 0, E s k(A v ) — 



A 



AAoA^llogA,,! 1 / 3 , where A 



Q4/3 5/6 



25 

Jq^O T siL-iQL-iv | lug, | ' , miicic /i - 2 n/6 

Experiment. We now briefly address experimental 
manifestations of the charge 2e skyrmions. The most 
direct way to observe the skyrmions in by the STM, 
which allows one to study the properties of the individual 
skyrmions, as well as the skyrmion configuration at finite 
density. For an individual skyrmion, STM may be used 
to study the dependence of the skyrmion size (|2H)) on the 
valley Zeeman interaction; the latter can be tuned by 
gates. Furthermore, in the vicinity of even filling factors, 
\v — 2M | <g; 1, there is a finite density of skyrmions in the 
system. When the skyrmion density n is low, such that 
the distance between skyrmions is much larger than the 
skyrmion size (|20[) . n" 1 / 2 3> a, the skyrmions form a tri- 
angular Wigner crystal. At larger densities, ttT 1 / 2 ~ a, 
when the skyrmions start to overlap, we find [23| that, 
similarly to the spin QHFM case [24| , the skyrmions rear- 
range into a square lattice. Such behavior will result in a 
periodic modulation of the local density of states, which 
can be measured by an STM and used to determine the 
skyrmion lattice symmetry. 
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